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Abstract. Context. Warps occurring in galactic discs have been studied extensively in HI and in the optical, but rarely in the 
near-infrared (NIR) bands that trace the older stellar populations. 
Aims. We provide NIR data of nearby edge-on galaxies, combined with optical observations, for direct comparison of the 
properties of galactic warps as a function of wavelength, and calculate warp curves for each galaxy and obtain the characteristic 
warp parameters. We discuss these properties as possible constraints to the different mechanisms that have been proposed for 
the development and persistence of galactic warps. 

Methods. We observed 20 galaxies that were selected from a statistically complete diameter-limited subsample of edge-on disc 
galaxies. We used the Cerro Tololo Infrared Imager (CIRIM) at the CTIO 1.5m Ritchey-Chretien telescope to acquire the NIR 
data. We used the 1.54m Danish and 0.92m Dutch telescopes at the European Southern Observatory's La Silla site for our 
, optical observations. 

Results. Our results show that 13 of our 20 sample galaxies are warped, with the warp more pronounced in the optical than 
at NIR wavelengths. In the remaining seven galaxies, no warp is apparent within the limitations of our automated detection 
method. The transition between the unperturbed inner disc and the outer, warped region is rather abrupt. SO galaxies exhibit 
Q\ • very small or no warps. The magnetic model remains one of a number of interesting formation scenarios. 
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1. Introduction therefore pay special attention to SO galaxies. Reshetnikov et 

rrt ■ .... , ,. al. (2002) find that warps were also common in the past, with 

Most spiral galaxies, including our own, exhibit warped discs. ., , - 

, , , , . r , , , even greater magnitudes at z ~ 1 . 
This has long been known based on observations or both the 

, , t , /y tt\ wo • imr n Edge-on galaxies are usually chosen to study warps, al- 

extended neutral-gas (HI) component (Sancisi 1976; Bosma ° ° J J . 

lrio1 t, . ,„ nn „ , „ . , .„„„,, , . though dynamical warps in the HI distribution of galaxies of 

1981; Bnggs 1990; Garcia-Ruiz et al. 2002b) and optical b . r b 

starlight (Sanchez-Saavedra et al. 1990; Florido et al. 1991; "folate indination have also been reported e^g in M83 

Reshetnikov & Combes 1998; Ann & Park 2006). The first de- T^S?^ NG ?" ' ^ 



NGC7331 (Bosma 1978), M31 (Brinks & Burton 1984), and 
M33 (Corbelli et al. 1989). Garcfa-Ruiz et al. (2002b) analysed 
HI observations and find that all edge-on galaxies that have an 
extended HI disc (with respect to the optical component) are 
warped. 



tections of warped discs came from 21 cm line observations of 
our own Galaxy (Burke 1957; Kerr 1957; Burton 1988). Based 
on a sky survey covering the Northern Hemisphere, Sanchez- 
Saavedra et al. (1990) first reported the high frequency of op- 
tical warps in disc galaxies, which they later confirmed in 

Sanchez-Saavedra et al. (2003) based on an extended catalogue Ann & Park ( 2006 ) find that 73% of the 325 S alaxies in 

that also included the Southern Hemisphere. These studies took ^ well-defined sample of almost perfectly edge-on galax- 

into account that warps can only be detected when their orien- les < a /* > 9 ' 5 ' wher « a > h 18 the major-to-minor axis ratio at 

tation (projection) is favourable (i.e., not in the line of sight). ^ = 25 ma § ^sec" 2 ) exhibit optical warps. Recently, Reyle 

A noteworthy and important result, established in the later pa- et aL ( 2009 ) have found that a wal ? is P resent in th * stellar ' 

per, is that they did not find any warped SO galaxies. Here, we dust ' and 8 as discs of the Mi]k y Wa ^ which are a11 asymmenic 
and characterised by a similar line of nodes. The Milky Way's 
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the stellar warp is significantly smaller (by a factor of approx- 
imately two). They conclude that this comparison shows that 
the different components react differently to the forces respon- 
sible for the origin of the warp. The most important additional 
information that we provide here (not yet been discussed by 
other authors) is a multiwavelength study of galactic warps in 
one NIR and three optical filters. 

Considerable effort has gone into establishing the physi- 
cal origin and stability of warps. Numerous mechanisms have 
been proposed and explored, but they have thus far not pro- 
vided a definitive and satisfactory physical explanation (see, 
e.g., the reviews by Binney 1992; Battaner et al. 1997; van der 
Kruit 2007). Kerr (1957), Hunter & Toomre (1969), Weinberg 
(1998), Weinberg & Blitz (2006), and others have suggested 
that warps can result from tidal interactions with a satellite or 
neighbouring galaxy. This hypothesis can explain warps nei- 
ther in isolated galaxies (Sancisi 1976; Tubbs & Sanders 1979; 
Sparke 1984) nor in the Milky Way, if the latter was pro- 
duced by an interaction with the Large Magellanic Cloud (see 
Garcfa-Ruiz et al. 2002a). Hunter & Toomre (1969) and Dekel 
& Shlosman (1983) proposed that the Galactic warp is a re- 
sult of oscillations in the disc, either triggered by an interac- 
tion with another galaxy in the past or caused by the disc be- 
ing embedded in a nonspherical halo of dark matter. Sparke & 
Casertano (1988) showed that some discrete oscillation modes 
can survive. This theory is not without problems, however; for 
instance, Binney (1991) showed that the halo should then re- 
spond to its misalignment with the disc, thus destroying the 
warp in a few orbital periods. Alternatively, a nonlinear cou- 
pling between a spiral wave and two warp waves could create 
a warp (Masset & Tagger 1997). The efficiency of this mecha- 
nism is too low in the stellar disc except at its outer edge (near 
a Lindblad resonance), where the spiral wave slows down and 
is coupled efficiently to the warp waves. At those radii, the en- 
ergy of the spiral arm can almost be completely converted into 
transmitted and reflected warp waves, which can be observed 
as corrugations in the disc. 

Warping could thus be a natural response of the outer disc 
to a series of stimulations, and it seems that the responsible 
mechanism does not necessarily have to be the same in each 
galaxy. In any case, observations of large edge-on galaxies 
could contribute to the definitive determination of the dominant 
scenario. The intergalactic medium could also play an impor- 
tant role by accretion of matter as galaxies move along their 
paths (Kahn & Woltjer 1959; Lopez-Corredoira et al. 2002; 
Sanchez-Salcedo 2006). Battaner et al. (1990) explained warps 
by intergalactic magnetic fields, which would produce a di- 
rect distortion in the distribution of the gas. The stellar system 
would also indirectly reach a warped distribution because stars 
are formed from gas, but this would result in a difference be- 
tween young and old stars. Colour gradients in warps can thus 
provide an important clue. Based on this model, gas-poor SO 
galaxies should not exhibit warps. 

In this paper, we present new, deep NIR and optical obser- 
vations of 20 edge-on galaxies. NIR wavelengths have rarely 
been used to date to study warps. They are, however, essential 
for better tracing the stellar mass distribution, because at in- 
frared (IR) wavelengths (1 < X < 5/im) dust extinction is min- 



imised (the absorption, in magnitudes, in the K filter is 10% of 
that in the V band, A K /A V ~ 0.1). NIR emission is also a better 
tracer of old stellar populations, one of our major objectives in 
this study. Therefore, these wavelengths provide essential com- 
plementary information for studying stellar warps. 

2. Sample 

Our sample was selected by de Grijs (1997, 1998) from the 
Surface Photometry Catalogue of the ESO-Uppsala Galaxies 
(ESO-LV; Lauberts & Valentijn 1989). This catalogue contains 
large numbers of galaxies that have been selected and param- 
eterised uniformly. The galaxies were originally selected on 
ESO-Schmidt survey plates covering declinations S < -17.5°, 
excluding the area within 15° of the Galactic Equator, at 
Galactic latitudes \b\ < 15°. We selected 20 candidate galax- 
ies from a statistically complete diameter-limited subsample of 
edge-on disc galaxies (de Grijs 1997, 1998), with the following 
characteristics (see Table 1): 

- £>25(-B) > 2.2', i.e., with optical diameters (at a surface 
brightness level of fig = 25 mag arcsec~ 2 ) larger than 2.2'; 

- Isolated galaxies, i.e., without significant neighbours within 
a distance greater than or equal to five times their optical 
diameters. They should be classified as noninteracting and 
undisturbed; 

- Inclination, ; > 87°; 

- Morphological type range from SO to Sd, i.e., revised 
Hubble types later than T = -2 (spiral and lenticular galax- 
ies). 

The inclinations were determined assuming an intrinsic 
flattening (b/a) = 0.11 (de Grijs 1997, 1998), so that our sam- 
ple contains highly inclined galaxies in which warps should be 
easy to detect. This could introduce a bias against early types 
or galaxies with prominent bulges that should be taken into ac- 
count in the discussion of the warp phenomena. Using this ap- 
proach, we try to avoid the common confusion between spiral 
arms, corrugations, and warps. 

3. Observations and data reduction 

The sample was observed at 

- Near-infrared wavelengths, with the Cerro Tololo Infrared 
IMager (CIRIM) at the 1.5 m Ritchey-Chretien Telescope 
of the Cerro Tololo Inter-American Observatory (CTIO) 
on October 22- -27, 1998 (Prop0020). CIRIM is equipped 
with a 256x256-pixel HgCdTe NICMOS3 array (Rockwell 
Science) with a pixel scale of 1 .16". We observed our sam- 
ple galaxies through the K s h 0It filter (=. K s ), characterised 
by an effective wavelength A e g = 2.15//m (Wainscoat & 
Cowie 1992). We opted for this filter instead of the longer- 
wavelength K filter (/t e ff = 2.20/im) because the sky bright- 
ness in K s is lower than in the K band (for instance, at 
UKIRT, Hawaii, the sky brightness is 13.5 and 13.0 mag 
arcsec~ 2 in K s and K, respectively 1 ) and K s hon is almost 

1 ht tp : // www . j ach . hawai i . edu/UKIRT/instrument s/ui st/ 
imaging/imaging .html 
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Table 1. Global galaxy parameters from the ESO-LV Catalogue (Lauberts & Valentijn 1989). 
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49.0 
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-31:32:09 


3.2+1 


6.0 


141.3 


14.68±0.25 


+ 1.38 


4992 


31.0 


ESO460-G31 
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11 Heliocentric radial velocities from Mathewson et al. (1992). 



as little affected by dust as the K band, A Ki ~ 0.1 1 mag 
ait-mass -1 ; and 

- Optical wavelengths. We obtained our observations us- 
ing two telescopes of the European Southern Observatory 
(ESO) at La Silla, Chile, including the 1.54 m Danish tele- 
scope (equipped with a 108 lx 1040-pixel TEK CCD with 
0.36" pixel 1 ) and the 0.92 m Dutch telescope (equipped 
with a 5 12x5 12-pixel TEK CCD with 0.44" pixel 1 ), both 
in the standard Johnson B and V and Thuan & Gunn (1976) 
i filters. 

The observations and optical data-reduction techniques on 
which the results presented in this paper are based were de- 
scribed in detail in de Grijs (1998). Here, we focus on the re- 
duction of the IR data. 

To study stellar warps in disc galaxies, we need to obtain 
reliable data at the very low surface brightnesses in their out- 
ermost regions. Our observational data reached surface bright- 
nesses fainter than 26 mag arcsec 2 in the B band. This is only 
6% of the sky brightness (jUB,sky = 23 mag arcsec -2 ), so that 
we had to be particularly careful in measuring the sky bright- 
ness. The sky background in NIR passbands is much higher 
than in the optical (pK,sky = 13 mag arcsec -2 ) and can change 
significantly on temporal and spatial scales the order of two 
minutes and two arcmin, respectively. It was therefore neces- 
sary to obtain as many sky frames as object frames to perform 
sky subtraction with sufficiently high accuracy. We took sky 
and object frames alternately, using the same exposure time for 
each individual frame (60 s) and employing interframe shifts of 
~ 2 arcmin. An additional reason we needed to obtain separate 
sky frames is that the IR arrays used were not large enough 
to properly sample the sky on the science frames. Incorrect 



background subtraction can cause the sky level to be either 
underestimated (so that sky flux will contribute to the surface 
brightness at large galactic radii) or overestimated (resulting in 
a fictitious radial cut-off). Figure 1 illustrates the importance 
of careful background subtraction in the K s band, where the 
sky contribution is greater than in optical passbands. The fig- 
ure shows the major- and the minor-axis (vertical, z) surface 
brightness profiles of four galaxies, binned radially. The effects 
of oversubtraction can be seen clearly in the minor-axis surface 
brightness profiles: they show artificial cut-offs, and the nega- 
tive background values result in undefined surface brightnesses 
at these z heights. 

We used a special software package for our NIR data re- 
duction, customised to reduce and analyse very deep data of 
extended objects observed in the IR, as a set of iraf 2 tasks. 
The program removes spurious stars in sky frames, subtracts 
cleaned, scaled sky frames from object frames taken immedi- 
ately prior to or after the relevant object frame, and flatfields 
the differences. It uses a median sky frame to detect and re- 
move stars in the sky frames. 

We observed several standard stars from the 
SAAO/ESO/ISO Faint Standard Stars Catalogue (Carter 
& Meadows 1995) to calibrate our targets and used the 
Wainscoat & Cowie (1992) correction to convert K s to K 
magnitudes. These stars were observed at least three times 
per night, at different airmasses, to calculate the prevailing 
atmospheric extinction with sufficient accuracy. We observed a 

2 iraf is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the US National Science Foundation. 
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Fig. 1. K s -band luminosity profiles of four galaxies. The lines in each panel represent, from top to bottom, -2cr (blue), -lcr (red), profile after 
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total of 21 standard stars. These observations reached a surface 
brightness of approximately 23 mag arcsec~ 2 in K^. 

4. Data analysis 

First, we aligned the images in all of the B,V,I, and K s filters. 
We took special care to subtract foreground stars. This is nec- 
essary and important because foreground stars projected close 
to a galaxy can contaminate the galaxy's luminosity distribu- 
tion and lead to erroneous results. We used a master mask con- 
structed by adding the individual masks for each of the pass- 
bands. In each mask, the stars were marked with circles of ra- 
dius around twice the FWHM to avoid contributions from any 
residual starlight. In some galaxies observed near the Galactic 
plane, this meant that a large number of pixels had to be 
masked. 

To determine the warp curve (tracing the position of the 
intensity maxima as a function of radius), we employed three 
different methods. 

1 . First, we calculated the intensity maxima using vertical cuts 
(perpendicular to the galactic plane). This resulted in very 
noisy profiles for the optical data and also sometimes for 
the NIR data, because the presence of a dust lane frequently 
resulted in a double-peaked luminosity distribution. 

2. Second, we calculated the luminosity distribution's first 
moment, with which we obtained reliable results (although 
it also depended on the position of the dust lane). This 
method was applied mainly at blue wavelengths to com- 
pensate for background stars. 

3. Third, we fitted Gaussian profiles. We fitted each vertical 
trace (parallel to the minor axis of the galaxy) in the se- 
lected zone to a Gaussian function. We considered this the 
best approach, because the contribution of the dust lane is 
less important in optical data than when using the other two 
methods. Gaussians can be matched to vertical profiles that 
contain a dip caused by a dust lane without problems in 
almost perfectly edge-on galaxies. 

4.1. Warp curves from the Gaussian fit 

As a first step, we rotated our galaxies (to align their major 
axes with the horizontal axis) using an iterative method until 
the position angle (PA) had been determined to an accuracy of 
approximately 0.05 degrees. We used the PA from the ESO-LV 
Catalogue as a starting point; we then iteratively fitted a straight 
line to the central region of the galaxy. Our results are shown 
in Fig. 2 as isophote maps of our 20 galaxies. 

To calculate the warp curves, we fitted Gaussians to the 
vertical profiles of the galaxies in each filter. We only con- 
sidered data with a signal-to-noise ratio greater than 3 (cr = 



Csys + cftd' wnere °"sys is the systematic error associated with 
the random background noise and cr std the standard deviation), 
and with the FWHM of the peak smaller than 100 pixels. The 
results of our analysis are presented in Fig. 3, where the warp 
curves are shown only for those values with error bars smaller 
than 0.5". This estimated error bar can be computed by scaling 
the standard deviation (lcr error) by the measured^ 2 value. 



4.2. Warp parameters 

It is convenient to use a geometric definition of a warp, e.g., for 
statistical studies. For an edge-on galaxy, the suitable coordi- 
nates correspond to the two directions contained in the plane of 
the sky: the direction defined by the major axis of the galaxy, 
x, and the direction of the rotation axis of the galaxy, y. The 
centre of the galaxy is the origin of coordinates. A warp curve 
is defined as the locus of points (x,,y,) tracing the distribution 
of the highest intensities, where y, is the position of the centre 
of the best-fitting Gaussian at a given x,. 

We derived the parameter co (defined in Jimenez-Vicente et 
al. 1997), 



CO 



1 r+L/2 
VJ-LI2 ' 



xydx, 



(1) 



where L is the full size of the edge-on galaxy (the diameter at 
3cr above the sky background). The absolute value of w is a 
measurement of the strength of a warp, while its sign distin- 
guishes between N-like (to > 0) and S-like (cj < 0) warps. 
These types are defined from the similarity of the warp shape 
and the type letters (see Fig. 4 for definition of types of warps). 
This definition is nondimensional and, therefore, independent 
of a galaxy's distance and size. It is also independent of the 
units used and does not strongly depend on the angular reso- 
lution of the observations (although a higher resolution would 
enable a more precise evaluation). For numerical purposes, we 
have 



1 



L/2 
i=-L/2 



(2) 



We can also define this parameter for both sides of our 
galaxies: w R for the right-hand side (Vz. 3 Zfio anc * Wl ^ or 
the left-hand side ('/l 3 ZjL-L/2 ^d')- Because of the flattening 
of the radial profile in the inner regions of our sample galax- 
ies, it is not critical to have a precise estimate of the position 
of the galaxy centre. Warps can appear in three flavours (cf. 
Lopez-Corredoira et al. 2008). The N- or S-like warps have 
o>r and o>l of the same sign, the signs of wr and <dl for U- 
like warps differ, and L-like warps apply to galaxies in which 
only one side is warped. The difference between N- and S-like 
warps is unimportant from a physical point of view. However, 
we have retained this information because it is needed, for in- 
stance, when considering the orientation of warps in a cluster 
of galaxies (see, e.g., Battaner et al. 1991). 

With this parameter we can, for example, objectively com- 
pare the appearance of a galaxy at several wavelengths to deter- 
mine whether a colour gradient exists within the warp. We can 
also detect warps that would otherwise not have been found. 
Warps are not always completely symmetric, in which case this 
parameter would hide the information of the clearly warped 
side of the galaxy. To avoid this problem, we calculated the 
warp parameters on either side of our galaxies independently, 
along with a s , the degree of asymmetry (Sanchez-Saavedra et 
al. 2003), 



I^r - fad 



(3) 
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S-shaped warp L-shaped warp 



Fig. 4. Definitions of angles, parameters, and types of warps. 



to the mean plane, and ends in the opposite hemisphere. The 
Milky Way's warp in the Southern Hemisphere is an exam- 
ple of this type of behaviour (Porcel et al. 1997). Florido et 
al. (1991) also report such 'elbow' -type warps. This effect may 
come from the presence of a more strongly warped dust lane 
(than the stellar disc) or be an intrinsic effect associated with 
warps. For such galaxies, a four-parameter fit would have pro- 
vided a better description. However, since an excessive number 
of fit parameters renders any interpretation of simple warps un- 
clear, we prefer to use three-parameter fits. 



One of the most important errors in estimating to arises 
from imperfectly determining the PA of the major axis, result- 
ing in an error 9 in the PA. Since 6 + 0, we introduce the error 



Aw 



= hL x{x 



i , e 

tan6»)dx « — T L 6 « -. 

3L 3 3 



(4) 



As 6 is, in general, approximately 0.05° (9 x 10~ 4 radians), we 
estimate Aw to be of the order of 3 x 10~ 4 rad. A large warp can 
have a value of 10 x 10~ 4 and a barely perceptible warp 5 x 1CT 4 
rad. In the remainder of this paper, we consider a galaxy as 
warped if its calculated cj is clearly greater than the associated 
error. (See Table 2 for our results.) 

For those galaxies showing an appreciably warped disc, we 
fitted a number of additional parameters to the warp curves. 
We tried to use the function proposed by Jimenez- Vicente et al. 
(1997), but it had too many fitting parameters so the fits were 
resticted to the simple description: 



\C(jc- 



\x\ < \A\ 
A) \x\ > \A\ 



(5) 



This function reproduces the warp's shape as a first approx- 
imation; i.e., it is flat up to a point and then deviates from the 
symmetry plane until it asymptotically reaches a new direction. 
The interpretation of the parameters A and C is as follows. 

- A is the starting point of the warp. It has dimensions of 
length. 

- C is the (nondimensional) value of the asymptotic slope. 

Throughout the rest of the paper, we refer to parameters 
A and C as those obtained by fitting the experimental points 
to Eq. 5. These values are shown in Table 3 for the differ- 
ent bands. In addition, we use the angles a and /3, which are 
most widely used to represent the warp amplitude (Sanchez- 
Saavedra et al. 2003) Here, a is the angle between the line con- 
necting the galaxy's centre and the outermost detected point, 
and the disc's major axis, while ft is the angle between the disc's 
major axis and the line from the outermost detected point to the 
point where the warp starts (A): see Fig. 4. 

We obtained the same result for the warp parameter w, 
whether it was calculated starting from the galaxy centre or 
from the starting point of the warp, A. This confirms that it 
is not critical to have a precise estimate of the position of the 
galaxy centre because of the radial flattening in the inner re- 
gions of our galaxies. In some galaxies, such as ES0142-G24 
or ES0157-G18, the warp sets off in one direction, turns back 



5. Analysis 

We consider a galaxy as warped if it is warped (according to 
the definition in Sect. 4.2) in any of the four observed bands. 
Based on this assumption, we find that 13 of our 20 sample 
galaxies (65%) are warped. This fraction supports previous re- 
sults based on larger galaxy samples with similar optical char- 
acteristics to our sample (Sanchez-Saavedra et al. 1990, 2003; 
Reshetnikov & Combes 1999). Given the difficulties of observ- 
ing warps with their line of nodes perpendicular to the line of 
sight, this high fraction of warped discs indicates that essen- 
tially all spirals are warped. Optical warps are less percepti- 
ble than HI warps, but at both wavelengths the frequency of 
warps is the same, i.e., close to 100% (Sancisi 1976; Garcfa- 
Ruiz 2001; Garcfa-Ruiz et al. 2002b; van der Kruit 2007). 

As expected, the N- and S-like warp frequencies are sim- 
ilar, since these characteristics depend on the observer rather 
than the galaxy. We observe six galaxies with N-like and four 
galaxies with S-like warps. The frequency of asymmetric warps 
is lower: ES0416-G25 exhibits a U-like warp, while ESO240- 
Gl 1 and ES0555-G36 show L-like warps. Therefore, the most 
frequent warp morphologies are N- and S-like warps (see also 
Ann & Park 2006). 

The distribution of warped galaxies as a function of mor- 
phological type is shown in Fig. 5, where the hatched area 
shows all galaxies and the cross-hatched area covers warped 
galaxies. The warp frequency does not seem to depend signifi- 
cantly on galaxy morphology, except for SO and very-late-type 
galaxies, which both yield lower frequencies. 

A noticeable result from Sanchez-Saavedra et al. (2003) 
was that SO galaxies are generally not warped. This places 
strong constraints on any hypothesis proposed to explain 
warps. Among our sample galaxies, ES0311-G12 has a small 
warp in both the optical and the NIR, but both of the other SO 
sample galaxies, ES0358-G26 and ES0358-G29, do not ex- 
hibit any detectable warps. Although our sample only contains 
three SO galaxies, so that we cannot draw any statistical con- 
clusions on that basis, our results are not in contradiction to 
previously published results. 

Since we have data in both optical and NIR bands, we 
can compare both, enabling us to either detect or reject colour 
dependences in the warps. This comparison is difficult, how- 
ever, because of inherent uncertainties and because the dif- 
ferences as a function of wavelength are not very large. Five 
of our galaxies (ES0235-G53, ES0288-G25, ES0311-G12, 
ES0416-G25, and ES0555-G36) show similar warp curves in 
all four observed bands, one galaxy (ES0531-G22) is warped 
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Table 2. Warp parameters for our four filters (xlO 4 ) and degree of asymmetry for each warped galaxy, with warp values listed in boldface. 



Galaxy 


B filter 


V filter 


/ filter 


K s filter 


oj 








to 








to 








to 








ESO026-G06 


1±3 


1±3 


-0±3 




0±2 


1±2 


-1±2 




2±3 


2±3 


1±3 




-3±3 


-2±3 


-1±3 




ESO033-G22 


-1±4 


0±4 


-1±4 




-3±6 


-1±6 


-2±6 




-1±4 


-1±4 


0±4 




3±4 


1±4 


1+4 




ES0142-G24 


3±2 


4±2 


-1±2 


1.0 


3±2 


4±2 


0±2 


1.0 


4±2 


4±2 


0±2 


1.0 


2±2 


1±2 


0±2 


- 


ES0157-G18 


-4±3 


-2±3 


-1±3 


- 


-4±2 


-3±2 


-1±2 


-1.0 


0±3 


-1±3 


1±3 


- 


1±3 


1±3 


0±3 


- 


ESO201-G22 


4±2 


3±2 


1±2 


1.0 


2±2 


3±2 


-1±2 


1.0 


2±2 


3±2 


-1±2 


1.0 


2±3 


2±3 


1±3 




hSUZOZ-G35 


17±6 


12±6 


5±6 


1.0 


16±3 


12±3 


5±3 


0.4 


11±6 


9±6 


3±6 


1.0 


4±4 


4±4 


2±4 




ES0235-G53 


29±6 


10±6 


19±6 


0.3 


39±7 


16±7 


22±7 


0.2 


28±3 


11±3 


17±3 


0.2 


25±2 


10±2 


15±2 


0.2 


ESO240-G11 


-4±3 


-4±3 


0±3 


-1.0 


-2±1 


-3+1 


1 + 1 


-1.0 


0±2 


-1±2 


0±2 




0±2 


0±2 


0±2 




ES0288-G25 


-6±5 


1±5 


-8±5 


-1.0 


-7±2 


-2±2 


-5±2 


-0.4 


-6±3 


-1±3 


-5±3 


-1.0 


-5±1 


-3±1 


-2±1 


-0.2 


ES0311-G12 


5±1 


3±1 


1±1 


0.5 


5±1 


3±1 


2±1 


0.2 


S±l 


3±1 


2±1 


0.2 


4+1 


2±1 


1 + 1 


0.3 


ESO340-G08 


-1 + 1 


1±1 


-2+1 


-1.0 


-1±1 


1±1 


-2+1 


-1.0 


-1±2 


1±2 


-2±2 




-1±1 


-1±1 


-1±1 


0.0 


ESO340-G09 


-2±7 


-6±7 


4±7 




4±5 


-1±5 


5+5 




1±8 


-6±8 


7±8 




9±9 


2±9 


8±9 




ES0358-G26 


-1 + 1 


-1±1 


-1 + 1 




0±1 


0±1 


0±1 




0±0 


0±0 


0±0 




2±2 


1±2 


2±2 




ES0358-G29 


0±1 


0±1 


0±1 




1±1 


1±1 


1±1 




-1±2 


0±2 


0±2 




-2±2 


-1±2 


-1±2 




ES0416-G25 


0±1 


2±1 


-1 + 1 


1.0 


0±1 


2±1 


-1 + 1 


1.0 


-2±2 


1±2 


-4±2 


-1.0 


-4±2 


-1±2 


-3±2 


-1.0 


ESO460-G31 


8±9 


1±9 


7 ±9 




4±4 


2±4 


3±4 




7±8 


1±8 


6±8 




0±3 


-3±3 


3±3 




ESO487-G02 


-1±3 


0±3 


0±3 




-1±2 


0±2 


-1±2 




1±2 


0±2 


1±2 




4±4 


0±4 


4±4 




ES0531-G22 


-20±6 


-14±6 


-6±6 


-0.4 


-16±4 


11±4 


-5±4 


-0.4 


11±4 


-4±4 


-7±4 


-0.3 


-6±2 


-2±2 


-4±2 


-0.3 


ES0555-G36 


-10±4 


0±4 


11±4 


-1.0 


-12±8 


-1±8 


-11±8 


-1.0 


-12±5 


-1+5 


11±5 


-1.0 


-10±5 


-2±5 


-9±5 


-1.0 


ES0564-G27 


3±3 


0±3 


2±3 




4±3 


2±3 


3±3 


1.0 


3±2 


1±2 


2±2 


1.0 


-5±1 


-3+1 


-2±1 


-0.2 




-2 2 4 6 8 



Morphological type 

Fig. 5. Morphological-type distribution of our warped galaxies. 
Hatched area: All galaxies. Cross-hatched area: Warped galaxies. 

in all four filters but exhibits a clearly weaker IR warp, one 
galaxy (ES0564-G27) is warped in the V,I, and K s filters, 
and six galaxies (ES0142-G24, ES0157-G18, ESO201-G22, 
ESO202-G35, ESO240-G1 1, and ESO340-G08) are warped in 
the optical bands but not in the IR. 

Since NIR profiles reach, in general, smaller radii than their 
optical counterparts, it is not straightforward to compare the 
warp distortion at different wavelengths. To do so, we com- 
pared the value of the fitted warp curve's y at the same radial 
distance in each band, for which we adopted the last measured 
point in 7£ s (i.e., the h parameter of Ann & Park 2006): see 
Table 3, column (3). Based on this table, we appreciate that the 
optical B and V bands do not exhibit significant differences. 




-4 -2 2 4 6 -2 2 4 6 

ly(B)l - ly(Ks)l ly(V)l - ly(Ks)l 



Fig. 6. Histogram of y(B) -y(K s ) and y(V) -y(K s ) at the last measured 
point in K s . 

However, we find a difference between V and K s : y(V) is usu- 
ally greater than y(X s ) (in 14 of 22 cases; see Fig. 6). On the 
other hand, y(K s ) is greater than y(V) for only six of the 22 
cases. This difference can be quite large. The most extreme case 
is ESO240-G1 1 (right-hand side), where the y(V)-z(K s ) = 5.5 
". We find a marginal tendency for NIR warps to have lower 
absolute values of the fit parameter C. Our attempts at statisti- 
cally detecting a colour dependence of the observed warps are, 
therefore, inconclusive. A better understanding would require 
more studies of this type. If the colour dependence were even- 
tually confirmed, mechanisms acting directly on the gas as the 
perturbing agent should be favoured. A colour difference has 
already been found by Florido et al. (1991). They also noticed 
that the dust lane in NGC 4013 is much more strongly warped 
than its optical disc. Recently, a study of 2MASS star counts in 
the Milky Way (Reyle et al. 2009) confirmed that the Galactic 
warp is less obvious in stars than in the gas and that the disc's 
shape is different at negative and positive longitudes. However, 
other authors have not been able to confirm such a difference 
(e.g., Momany et al. 2006 for the Milky Way). 

Figure 7 shows the behaviour of the warp angle Basa func- 
tion of passband. This angle represents the strength of the warp, 
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Table 3. Warp parameters as a function of passband. 



Galaxy name 






v (±0.5) 


A(±l) 


C(±0.01) 




Side 


ria^i ±1 (") 

(3) 


B 

(4) 


V 
(5) 


I 

(6) 


(7) 


B 

(8) 


V 
(9) 


I 
(10) 


(11) 


B 

(12) 


V 
(13) 


I 

(14) 


(15) 


ES0142-G24 


Right 


1 10 


-1 


-1.5 


1 





108 


106 


103 


- 


0.17 


0.12 


0.27 


- 




Left 


103 





0.5 








-105 


-103 


-104 


- 


0.10 


0.04 


0.08 


- 


ES0157-G18 


Right 


90 


-1.5 


-2 


-1 





70 


68 


- 


- 


-0.06 


-0.07 


- 


- 




Left 


92 


0.5 


0.5 





1 


-81 


-115 


- 


- 


-0.05 


0.04 


- 


- 


ESO201-G22 


Right 


71 





0.5 


0.5 





46 


63 


29 


- 


0.07 


0.10 


- 


- 




Left 


72 


1 


1 


0.5 





-66 


-86 


-63 


- 


0.05 


0.21 


0.02 


- 


ESO202-G35 


Right 


82 


2 


2 


2 





30 


33 


43 


- 


0.90 


0.07 


0.05 


- 




Left 


72 


-2 


-2 


-2 


-1.5 


-41 


-44 


-45 


- 


0.07 


0.07 


0.09 


- 


ES0235-G53 


Right 


70 


7 


6 


5 


6 


27 


36 


30 


30 


1.52 


0.16 


0.15 


0.39 




Left 


70 


-7 


-5 


-6 


-4 


-11 


-35 


-35 


-39 


0.11 


0.16 


0.1 1 


0.10 


ESO240-G1 1 


Right 


180 


-5.5 


-5.5 


-3 





101 


132 


- 


- 


- 


- 


- 


- 




Left 


165 


-0.5 


-1.5 


-1 







-103 






0.00 


0.02 






pcmofi no's 


Ri a ht 


o / 


-0.5 


-1 


-3.5 


-2.5 


65 


60 


60 


43 


-0.08 


-0.09 


-0.49 


-0.04 




Left 


73 


2.5 


2 


1 





-26 


-28 


-29 




-0.03 


-0.05 


-0.02 




ES0311-G12 


Right 


118 


1.5 


1 


1 





46 


32 


47 


32 


0.03 


0.02 


0.02 


0.02 




Left 


102 











-49 


-75 


-63 


-60 


0.01 


0.06 


0.03 


0.02 


ESO340-G08 


Right 


37 














64 


59 






-0.08 


-0.12 








Left 


39 














-50 


-50 






-0.05 


-0.05 






ES0416-G25 


Right 


63 


1.5 


1.5 


1.5 


-3 


51 


51 


51 


59 


1.00 


0.68 


0.61 


-0.68 




Left 


69 










1 






-46 


-54 






-0.04 


-0.06 


ES0531-G22 


Right 


88 


-4 


-3.5 




-2.5 


58 


44 


58 


44 


-0.14 


-0.12 


-0.11 






Left 


69 


2 


1.5 




2 


-39 


-48 


-12 


-57 


-0.05 


-0.04 


-0.04 


-0.14 


ES0555-G36 


Right 


52 
















53 


50 








-0.07 






Left 


69 


2.5 


2 




1.5 


-28 


-25 


-28 


-21 


-0.04 


-0.02 


-0.04 


-0.02 


ES0564-G27 


Right 


117 





0.5 


0.5 


-1 








48 






-0.04 


-0.04 




Left 


83 





-1 


-0.5 


-0.5 




-25 


-70 


-37 




0.01 


0.02 


-0.02 



(3) Radius of the last measured point in K s (arcsec). (4), (5), (6) and (7) y at the same point (arcsec). 
(8), (9), (10) and (11) Warp starting points (arcsec). (12), (13), (14) and (15) Asymptotic slope. 
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Fig. 7. Distribution of warp angle as a function of revised Hubble type. 



i.e. how much the outer disc deviates from the plane defined 
by the inner (unwarped) disc. Its highest value is obtained for 
galaxies with revised Hubble type T = 3, corresponding to Sb. 
For these galaxy types, a tendency exists toward lower values 
of B in the NIR filter than in the B band. The starting point of 
the warp, A, does not appear to be a function of wavelength. 

Kregel et al. (2002) analysed and used the same optical ob- 
servations discussed here to study the flattening and truncation 



of the stellar discs. We take their measurements of the trunca- 
tion radius (see van der Kruit 2007 for definition), R nrdx , in the 
/ band for 1 1 of our warped galaxies. It would be very useful 
to compare this truncation radius with the starting radius of the 
warp (A), because the relation between warps and disc trunca- 
tions in edge-on galaxies (van der Kruit 2007), if any, is un- 
clear. Obviously we obtain from this comparison that the warp 
radius, A, is approximately equal to 0.8 R max . Van der Kruit 
(2007) concludes that when an HI warp is present, it starts at 
approximately l.l/? ma x- Therefore it appears that optical warps 
always start earlier than HI warps (HI warps start at approxi- 
mately 1 .4 A), although we note that the low resolution of the 
HI data makes it difficult to detect low-amplitude warps. 

Kregel et al. (2004) also presented HI synthesis observa- 
tions of 15 edge-on galaxies. They have seven galaxies in com- 
mon with our sample, of which ESO487-G02 was observed 
with insufficient sensitivity, and five galaxies are found to have 
a warped HI disc. We now analyse the objects in common with 
Kregel et al. (2004) for which warps have been detected either 
by them, by us, or by both. They find that the neutral hydrogen 
distribution in ES0142-G24 appears warped on both sides of 
the galaxy, following a similar pattern to the stellar disc. In our 
optical observations, the warp starts at 106 " in the south and 
104 " in the north, but on the basis of our NIR observations, 
we could not detect the same behaviour. Instead, we observe a 
corrugated disc. 

Kregel et al. (2004) find that ES0157-G18 is lop- 
sided, which may be related to the small companion, 
APMBGC157+052+052, although this galaxy had no previ- 
ously determined redshift. They did not detect any clear HI 
warp. We also detected some differences between both sides 
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of this galaxy. Our warp curves exhibit a similar behaviour in 
all optical bands, while in the NIR the warp is not clear. 

ESO201-G22 shows a similar warp curve in HI and in all 
optical bands, but it is uncertain in the NIR, where the data are 
not as deep as in the optical. The /T s -band warp curve finishes 
at a smaller radius than in the optical bands. The bright star on 
the eastern (left-hand) side of the galaxy hinders our analysis 
of this region. 

ESO240-G11 is very interesting because a region in the 
galaxy's northwest (right) is clearly warped in the B and V 
bands, but the warp is not as marked in the NIR. Kregel et 
al. (2004) detected warping of the HI layer and they point out 
that the HI warp can be traced by faint emission in the /-band 
image. In the optical bands, the dust lane is not completely co- 
incident with the galactic plane, which means that this galaxy 
is either not perfectly edge-on or that it has a warped or asym- 
metric bulge, which also results in a rather large error bar. On 
the southeastern (left-hand) side, the warp radius is 103 + 2 " 
with B = -2.3°, and on the northwestern side this radius is 
132.3 + 0.6 " and 8 - 63.7°, so this galaxy exhibits an L-like 
warp. The 'elbow' where the warp begins in a given direction 
turns back to the mean plane and ends in the opposite hemi- 
sphere should show the beginning of the spiral arms. A more 
detailed study of corrugations in the galaxy's stellar and dust 
discs would also be very interesting. 

The U-like type galaxy ES0416-G25 shows a similar warp 
curve in all bands, including in HI and the NIR. 

The final galaxy also studied by Kregel et al. (2004) is 
ES0564-G27. This is a clearly corrugated Sc galaxy with some 
bright foreground stars projected onto the major axis, so it 
is very difficult to obtain robust conclusions about the stellar 
warp. Kregel et al. (2004) find that the neutral hydrogen is 
rather symmetric and extends further than the stellar disc. It 
is warped on both sides, with the warp apparently starting be- 
yond the stellar disc. We cannot compare this with our results 
because of the large uncertainties. 

Then we can conclude from this comparative study that 5 
out of 6 galaxies have an HI warp (only ES0157-G18 does not 
show any clear HI warp), and these 5 galaxies have an optical 
warp, while in the NIR the presence of a warp is unclear or 
absent in all of them. This would indicate that the old star pop- 
ulation is less warped than the gas and the more recently born 
stars. 

We calculated a s - see Eq. (3) - for all warped galaxies, 
and the results are shown in the Table 2. The asymmetry in 
the warp angles seems related to the warp's intrinsic proper- 
ties. Garcfa-Ruiz et al. (2002b) show that the largest asymme- 
tries are observed in galaxies that display obvious tidal features. 
We can only show, based on our sample, that the three galax- 
ies that exhibit asymmetries, ES0235-G53, ES031 1-G12, and 
ES0531-G22, all show a large warp amplitude, which supports 
the positive correlation found by Garcfa-Ruiz et al. (2002b). 
However, none of the galaxies in our sample are involved in 
strong tidal interactions. 



6. Discussion and conclusions 

Are all warps produced by the same mechanism? For example, 
gravitationally driven warps should exhibit no differences be- 
tween gaseous and stellar morphologies, as opposed to mecha- 
nisms acting on the gas (e.g., magnetically or accretion-driven 
warps). Are the warp curves the same for different colours; 
i.e., is there any difference in shape and magnitude at differ- 
ent wavelengths? 

The present study suggests that, in some of the observed 
galaxies, the warp is smaller at NIR wavelengths than in the 
optical. This - combined with the important result of Sanchez - 
Saavedra et al. (2003) that lenticulars do not exhibit warps 
(which is also compatible with the results in this paper) - sug- 
gests that gas seems to be a necessary ingredient in the devel- 
opment of warps. We have not been able to robustly confirm 
the absence of warps in SO-type galaxies because our sample 
only contains three SOs. Therefore, we cannot draw any sta- 
tistical conclusions on that basis, but our results (at least) do 
not contradict previously published results that most (if not all) 
lenticulars are not warped. This must be explained by some 
theoretical model. Models based on gravity alone might have 
significant difficulties to explain these observations. A recent 
study using 2MASS star counts in the Milky Way (Reyle et 
al. 2009) shows that the warp is less pronounced in stars than 
in the gas, and it also shows different disc shapes at negative 
and positive longitudes. The magnetic model of Battaner et al. 
(1990) is consistent with all of these observations. 

Our results show that NIR warps exist with a frequency that 
is nearly the same as that of optical warps. This suggests that 
warps are long-lived structures that were formed in the early 
stages of the evolution of most galaxies. Since warps can be 
dissipated on a time scale close to the galaxy's rotation pe- 
riod, this indicates that the responsible mechanism, whatever 
it may be, seems to be acting permanently. The intergalactic 
magnetic-field model remains an interesting possibility. It is 
well-known that a warp will disappear on a time scale of ap- 
proximately 2 Gyr (Binney 1992). However, in the magnetic 
model of warps, the field configuration is not an initial con- 
dition but extragalactic fields act permanently. The extragalac- 
tic magnetic-field lines could be frozen in in the intergalactic 
medium, therefore corotating with galaxy clusters and present- 
ing a relatively stable configuration of the field with respect to 
the warped galaxy. 

The magnetic model was designed to explain N- and S-like 
warps (m = 1 modes), the most frequently occurring type (51% 
out of 73% of warped discs in the large sample of 325 edge-on 
galaxies of Ann & Park 2006). This model does not manage to 
explain m = warps (U-shaped profiles). They should be inter- 
preted as coming from gradients in the extragalactic field, with 
characteristic lengths close to a galaxy's size. Generations of 
asymmetric warps through other mechanisms have been con- 
sidered by Saha & Jog (2006) and Lopez-Corredoira et al. 
(2002). Other mechanisms have been proposed to induce warp- 
ing. In fact, the variety of warp morphologies could come from 
different causes that may act differently from galaxy to galaxy. 
However, observations of NIR warps and direct comparison 
with their optical counterparts could provide some constraints. 
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The frequency of warped discs in our sample is very high 
(65%). Given the difficulty of detecting warps with their lines 
of nodes in the plane of the sky, this suggests a fraction of 
nearly 100% for the frequency of warps in spirals. We have also 
been able to confirm the finding by Reshetnikov & Combes 
(1999) that warps are equally present in all types of spirals, 
T>0. 

The abrupt break between the undisturbed inner disc and 
the warped region (with a roughly constant slope) suggests 
that these two regions have different formation and evolu- 
tionary histories (e.g., van der Kruit 2007). However, van der 
Kruit concludes that HI warps start at l.lR max . If this were 
the case, optical and NIR warps could not be observed. That 
they are indeed observed allows us to conclude that R max > A. 
Nevertheless, we cannot discard a possible relationship be- 
tween truncations and warps. It appears that optical warps al- 
ways start closer in than HI warps, although we note that the 
low resolution of the HI data makes it difficult to detect low- 
amplitude warps. 
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Fig. 3. Left: Warp curves and right: Contour maps, for the different filters of our 20 galaxies. The warp curve (green line) is put on top only for 
those values with an error bar less than 0.5". The red line shows the fitted warp curve. The isophotes are equidistant (in units of n*3cr equiv. to 
a step of +0.75 mag/arcsec 2 ) starting at a level of about 3cr above the sky background. The galaxy is rotated, but N is closer to the top and E to 
the left. The scale on the y-axis is enlarged to show a better curve. 
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Fig. 3. continued 
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